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Abstract 


Within the one-loop approximation of baryon chiral perturbation theory we 
calculate all one-pion and two-pion exchange contributions to the nucleon- 
nucleon interaction. In fact we construct the elastic NN-scattering ampli¬ 
tude up to and including third order in small momenta. The phase shifts 
with orbital angular momentum L > 2 and the mixing angles with J > 2 
are given parameterfree and thus allow for a detailed test of chiral symmetry 
in the two-nucleon system. We find that for the D-waves the 27r-exchange 
corrections are too large as compared with empirical phase shifts, signaling 
the increasing importance of shorter range effects in lower partial waves. 
For higher partial waves, especially for G-waves, the model independent 
27r-exchange corrections bring the chiral prediction close to empirical NN 
phase shifts. We propose to use the chiral NN phase shifts with L > 3 as 
input in a future phase shift analysis. Furthermore, we compute the irre¬ 
ducible two-pion exchange NN-potentials in coordinate space. They turn 
out to be of van-der-Waals type, with exponential screening of two-pion 
mass range. 
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I. INTRODUCTION 


The force between two nucleons is one of the fundamental problems in nuclear physics. 
Experimentally, the two-nucleon force is mapped out in nucleon-nucleon scattering, which 
is a purely elastic process up to nucleon laboratory kinetic energies of Ti ab = 280 MeV, 
the iViV7r 0 -threshold. Since both the target and the projectile have spin-1/2, there is a 
rich spin-structure in NN-scattering. Besides total and differential cross sections many 
independent spin-observables can be measured with polarized targets and/or projectiles. 
At present a huge body of NN-scattering data exists and this data base will increase and 
improve in the near future, in particular with experiments to be performed at the proton 
cooler facility COSY in Jiilich. 

The theoretical interpretation of the elastic NN-scattering data is done in terms of a 
phase shift analysis. The pertinent phase shifts are labeled by the spectroscopic notation 
2S+1 Lj, with S = 0,1 the total spin, J = 0,1,2... the total angular momentum and 
L = J — S, J, J + S the orbital angular momentum of the two-nucleon system. The total 
isospin of the NN-system / = 0,1 is readily determined by the Pauli exclusion principle, 
which demands / + S + L to be odd, in order to have an antisymmetric wave function in 
the combined isospin-spin-coordinate space. The spin-dependence of one-pion exchange 
(longest range NN-interaction) gives rise to a tensor interaction which causes a mixing of 
the (triplet) states with orbital angular momentum L = J — 1 and L — J + 1. In order to 
quantify this effect a mixing angle ej is introduced for each J in addition to the singlet 
and the three triplet phase shifts. The empirical mixing angles €j are, however, small and 
do not exceed 7 degrees up to the YY7r°-threshold. This weak mixing justifies the use of 
the spectroscopic notation 2S+1 Lj which presumes L to be a good quantum number. 

The single (weakly) bound state, the deuteron in the 3 5/-channel with total isospin 
1 = 0, shows up as a pole in the NN-scattering amplitude (on the physical Riemann-sheet 
of Ti a b ) slightly below threshold at Ti ab = —4.45 MeV, corresponding to minus twice the 
deuteron binding energy [lj. In the 1 S'o-channel with total isospin 1=1 the nucleon- 
nucleon potential is slightly too weak to form a bound state and a pole on the second 
Riemann-sheet results |2[], even much closer to threshold than the deuteron pole. Both 
the true bound state deuteron and the ’’anti-bound” state on the second Riemann-sheet 
slightly below threshold are responsible for the extraordinarily large S-wave NN-scattering 
lengths 0 (a( 3 S'i) = —5.42 fm, a( 1 S'o) = 23.75 fm) and the rapid change of the S-wave 
phase shifts. In all other partial waves the energy dependence is much more moderate. P- 
and D-wave phase shifts change by at most 30° from threshold up to the VATY-threshold 
at T iab = 280 MeV and F- and G-waves do not exceed values of 7° in this energy range Q. 
Thus with the exception of the two S-waves, the nucleon-nucleon interaction is actually 
rather weak in almost all channels. 

An accurate description of the NN-scattering data and deuteron properties has been 
found in terms of the one-boson-exchange model (OBE) and more refined the Bonn- 
potential |Q which includes multi-meson exchanges as well. The basic ingredients of the 
OBE model are the one-pion exchange, the exchange of the vector mesons p(770) and 
o;(782), and the exchange of a fictitious scalar isoscalar boson cr(550). The latter may be 
viewed as an effective description of correlated two-pion exchange in the isoscalar scalar 
channel. Such an interpretation of the cr(550) conies from the dispersion-theoretical stud¬ 
ies of the two-pion exchange [Q, which use the empirical information from the crossed 
7rV-scattering channel. The coupling constants of the vector and scalar mesons to the 
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nucleon are adjustable parameters. For each meson-nucleon vertex a ’’form factor” of 
monopole-type with adjustable cut-off in the range 1 GeV < A < 2 GeV is introduced. 
In the OBE model these form factors are needed to have convergent integrals in the 
Lippmann-Schwinger equation, which iterates the OBE potential to infinite orders. Even 
though the meson-nucleon form factors are often interpreted (in analogy with the electro¬ 
magnetic nucleon form factors) in terms of nucleon substructure, they do not represent 
truely physical quantities. First, they cannot be measured, and secondly they cannot even 
be defined in a model independent way. As off-shell matrix elements they depend on the 
actual form of the interpolating meson held, which is not unique for a composite particle. 
Putting aside the non-observability of meson-nucleon form factors, the OBE-model works 
very well. With about a dozen free parameters one obtains accurate fits to NN-scattering 
data and deuteron properties |5J. 

In recent years there has been activity to understand the NN-force in a more fun¬ 
damental way from the symmetries of QCD, in particular from chiral symmetry which 
governs low-energy strong interactions. In his seminal papers Weinberg 0 discussed the 
qualitative implications of chiral symmetry for the two- and multi-nucleon forces. Or¬ 
donez, Ray and van Kolck || applied the chiral Lagrangian to the two-nucleon problem. 
Using time-ordered perturbation theory and a Gaussian cut-off prescription to regularize 
loop divergences they calculated the two-pion exchange NN-potential. In addition more 
than 20 parameters related to short range NN-contact terms were introduced. The cut¬ 
off dependent potential was then transformed to coordinate space and used to solve the 
non-relativistic Schrodinger equation. A good fit was obtained for the deuteron properties 
and the phase shifts in the low partial waves (S, P, D) up to kinetic energies of Ti ab = 100 
MeV (partly also up to Ti ab = 300 MeV) ||. Of course in the presence of nearly 30 fit 
parameters the role of chiral symmetry in the NN-interaction is obscured. Furthermore, 
the use of a Gaussian cut-off is not consistent with the chiral power counting scheme 
underlying the whole approach 0. 

The purpose of this work is to explore the kinematical window in energy and angular 
momentum in which the NN-interaction is governed by chiral symmetry alone. Our 
physical motivation comes from the simple fact that the range of the force produced by 
the exchange of a hadronic state is inversely proportional to its mass. The longest range 
part comes from one-pion exchange and chiral symmetry manifests itself via the smallness 
of the pion mass, a feature which derives from the fact that the pion is an approximate 
Goldstone boson of spontaneous chiral symmetry breaking in QCD. The exchange of 
the next heavier hadronic state, consisting of two pions, is intimately related to ttN- 
scattering by crossing. In this reaction chiral symmetry leads to important dynamical 
constraints in the form of low-energy theorems [|9j. see also the recent work of [|10J, where 
many low-energy 7r N-o bservables have been investigated in chiral perturbation theory. It 
is now interesting to explore the structure of the 27r-exchange NN-potential when it is 
produced by a 7T ^-interaction with all chiral constraints built in. Naturally, one expects 
such an approach to work for the peripheral NN-partial waves where only the long and 
intermediate range components of the NN-force operate. 

In the present work we calculate the chiral two-pion exchange based on the same 
chiral nN -Lagrangian as used in 0. Our methods and aims are however different and 
complementary. First of all we will use covariant perturbation theory throughout and 
thus have to solve the problem of the so-called pinch singularity [|7j occuring in the planar 
box diagram. Secondly, we use throughout dimensional regularization, a prescription to 
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handle loop divergences which is in harmony with chiral symmetry and the chiral power 
counting scheme. Thirdly, we stay within the systematic expansion of chiral perturbation 
theory and calculate in one-loop approximation all contributions to the NN T-matrix up 
to and including third order in small external momenta. Consequently, we do not iterate 
to infinite orders. As a result we find that the phase shifts with L > 2 and mixing angles 
with J > 2 are given entirely parameterfree and thus allow for a rigorous test of chiral 
symmetry in the NN-interaction. 

The paper is organized as follows. In section 2, we briefly discuss the chiral nN- 
Lagrangian at next-to-leading order underlying our calculation. In section 3, we introduce 
the necessary formalism, the nucleon-nucleon T-matrix, the projection onto partial waves 
and the expressions for the perturbative phase-shifts and mixing angles. Section 4 presents 
closed form analytical results for the one-pion and two-pion exchange contributions to 
the NN T-matrix calculated in one-loop approximation. In section 5 the parameterfree 
results for the peripheral phase shifts (L > 2) and mixing angles (J > 2) are discussed. In 
section 6 we give explicit expressions for the irreducible two-pion exchange potentials in 
coordinate space and compare them with phenomenological ones. In section 7 we compare 
the isoscalar central amplitude at zero momentum transfer with OBE and finally section 
8 ends with a summary. In the appendix the construction of the anti-symmetric NN 
T-matrix satisfying the Pauli exclusion principle is given. 


II. EFFECTIVE CHIRAL LAGRANGIAN 


The tool to investigate the dynamical consequences of spontaneous and explicit chiral 
symmetry breaking in QCD is the effective chiral Lagrangian. It provides a non-linear 
realization of the chiral symmetry group SU(2)l x SU(2)r using the relevant low-energy 
effective degrees of freedom, the Goldstone pions and the nucleons. The effective chiral 
7T N - Lagrangian takes the following general form, 

£ = + • • • (i) 


where the superscript indicates the number of derivatives or small external momenta. 
In the heavy baryon formulation [|nj (which essentially corresponds to a non-relativistic 


treatment of the nucleons) the leading and next-to-leading order terms read |fl0| . p~2 


£S = JV(iD„- 


( 2 ) 


4 » = jv| 


D ■ D + i^-,{3 ■ D, u„] + 2d milt' + rf) 


2 M 


AM 


2 1 
+G - iuH + + H c * + 1bp • (a x 3) K 


8 M 


AM' 


( 3 ) 


Here, D M is the chiral covariant derivative acting on the iso-doublet 

nucleon field N and is an axial vector (matrix) quantity. The 577(2)- 

matrix U = £ 2 collects the Goldstone pion fields in the form U = l+ir-n/— 7f 2 /2/^ + ... 
with the pion decay constant = 92.4 MeV. The only other parameter occuring in 
is gA, the axial vector coupling constant of the nucleon. Since we are considering 
here exclusively the pion-nucleon vertex we employ the Goldberger-Trciman relation gA = 


A 











g-irNf-n/M together with the empirical value of g^N = 13.4 |H| and therefore use gA = 1-32 
throughout. M = 939 MeV is the (average) nucleon mass and m n = 138 MeV the 
(average) pion mass. 

The second order Lagrangian C^ N consists of two different types of terms. The first two 
terms in eq.(3) have coefficients 1/2 M and paI^M which are fixed by Lorentz invariance 
IT]. The remaining set of chiral invariant terms are accompanied by new low-energy 
constants ci, 02 , 03 , 04 . These four low-energy constants have been recently determined 
at one-loop order in ref. [jTDJ from a fit to nine (sub)threshold 7 riV-observables. Their 
values are C\ = —0.9 ± 0.1, c 2 = 3.3 ± 0.2, c 3 = —5.3 ± 0.2, C 4 = 3.6 ± 0.1, all given in 
GeV -1 . As shown in ref. [|TD| the A(1232)-resonance makes the dominant contribution to 
c 2 , c 3 and C4. The exchange of a scalar boson with mass and coupling constant similar to 
the fictitious cr(550) allows to explain c\ and also part of c 3 . Furthermore, C 4 receives a 
significant contribution from p(770) vector meson exchange. Note that the c ? ; are much 
larger than the typical l/2M-scale in eq.(3). Thus one expects the important 27r-exchange 
contributions to the NN T-matrix to be the ones proportional to q. This completes the 
description of the effective chiral 7 nY-Lagrangian. 


III. FORMALISM 

We are considering elastic nucleon-nucleon scattering N(p) + N(—p) —> N(p') + 
N(-p'), where ±p and ±p' are the initial and final nucleon momenta in the center-of- 
mass (CM) frame with p = \p\ = \p'\. The momentum transfer q — p' — p has the 
magnitude q — \q\ = 2p sin | = p^j2(1 — z). Here 9 is the CM scattering angle with 
z = cos 9. In terms of the nucleon laboratory kinetic energy T/ 0 &, the CM momentum p is 
given as p = sjT tab M/2. 

In the center-of-mass frame the elastic on-shell nucleon-nucleon T-matrix takes the 
following general form, 

T nn = Vc + Ti • T 2 Wc + [V5 + T\ ■ T2W5] CTi • cr 2 + [Vt + T\ • T2W7’] < 7 \ ■ q (J 2 • q + 

[Vso + Ti • t 2 Wso\ i(vi + 02 ) ■ (g X p) + [V Q + n ■ f 2 W Q ] d x ■ {q x p) a 2 ■ (q x p) (4) 

with the ten complex functions Vc(p, z),.. . , Wq(p, z ) depending on p and z. Here ofi and 
cx 2 are the spin-vectors of the two nucleons. The subscripts refer to the Central, Spin-spin, 
Tensor, Spin-Orbit and Quadratic spin-orbit components of the NN T-matrix. Each of 
these five components occurs in an isoscalar (V) and an isovector (W) version. Note that, 
for the purpose of having projection formulas as simple as possible (see eqs.(6-9)), we do 
not follow the usual conventions for the decomposition of T^n in all cases. In particular, 
we associate the tensor term with the helicity structure • qa 2 ■ q in our decomposition. 
The overall sign of T^n is chosen here as it follows directly from the Feynman rules for the 
S-matrix, with a positive T-matrix, positive phase-shifts and a positive scattering length 
corresponding to attraction. 

We remark that eq.(4) is not just a non-relativistic approximation, but results from 
the evaluation of the fully relativistic NN T-matrix in the center-of-mass system when 
reduced to two-component spinors. In principle the NN T-matrix consists of the total 
sum of all direct diagrams and those with crossed nucleon lines. In practice, there is 
however no need to evaluate the diagrams with crossed nucleon lines, since the Pauli 
exclusion principle gives a simple rule, namely L + S + I being odd, to select the allowed 
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NN-partial waves. Therefore we identify here Tnn in eq.(4) only with the sum of all 
direct diagrams. In the appendix it is shown how to construct from the sum of direct 
diagrams the complete anti-symmetrized form by application of Fierz-transformations and 
the substitution z —> — z. 

In order to compute phase shifts and mixing angles the matrix elements of Tnn hi the 
LSJ -basis are needed. To obtain these one first evaluates matrix elements in the helicity 
basis and then rotates with Wigner d-functions to the LSJ- basis. We have followed the 
detailed description in ref. ||14| 
with total isospin / = 0 , 1 , 


and find with the relevant linear combination in a state 


U K = Vk + (4/ - 3 )W K , (K = C, 5, T, SO, Q ) 


the following projection formulas: 

a) Singlet matrix element with S — 0 and L = J: 

(J0J\T nn \J 0J) = l f 1 dz\u c - 3U S - q 2 U T + p 4 (z 2 - 1 )U Q \Pj(z ). 

Z J-1 J 

b) Triplet matrix element with S — 1 and L = J: 

(J1J\T NN \J1J) = l - /;' dz{2p 2 [u so -U T + p 2 zUq ] (P J+1 (z) + Pj-i(z)) 

+ [U c + U S + 2p 2 {l + z)U T - Ap 2 zU so ~ p\3z 2 + 1 )U Q \ Pj(z)} . 

c) Triplet matrix elements with 5 = 1 and L — J ± 1: 

(J ± 1, 1J\T nn \J± 1,1 J ) = i £ dz{2p 2 [u so ± 2 T+i( Ut ~ P 2 zU q)\ P M 
+ [U C + U s +p 2 (p 2 { 1 - z 2 )U q - 2zU so ± ^-^{p 2 U Q - U T )) 

d) Triplet mixing matrix element with 5 = 1, L' = J — 1 and L = J + 1: 

V7Tip 2 r 1 


(5—1, IJ\T nn \J + 1,15) — 


dz{(U T -rrU Q )P J+1 (Z 


(5) 


( 6 ) 


(7) 


C;±i(*-)}. (8) 


(9) 

we 


\fj ( 2 J + 1 ) J -1 

+ [(25 - z{2J + 1 ))U t + P 2 zU q \Pj(z)} . 

Here, Pj(z) are ordinary Legendre polynomials of degree 5. In accordance with ref. |14[ 
did not use the more conventional tensor operator 3+i ■ q +2 • q — d\ ■ Jo in the decomposition 
of T nn , eq.(4), since it leads to unnecessarily complicated projection formulas. 

The phase shifts and mixing angles (in the standard convention of Stapp et al. |15| ) 
are then given perturbatively [|T 6 | 


as 


5lsj = ^ Re (LSJ\T nn \LSJ ) , ej = y-qy Re(5 - 1, 1J\T NN \J + 1,15) , (10) 

47 irj 


with the CM nucleon energy E = \JM 2 + p 2 . The kinematical prefactor comes from the 
proper relativistic flux and two-body phase space factors. Since we use chiral perturbation 
theory here, there is no exact but only a perturbative unitarity. Therefore the calculation 
of the phase shift according to eq.( 10 ) is valid only as long as the difference between 6 
and sin 5 cos 6 is small. This holds reasonably well for |<5| < 30°; thus with respect to this 
constraint all NN phase shifts and mixing angles, with the exception of the two S-wave 
phase shifts, can be calculated reliably in (chiral) perturbation theory. 


6 




























IV. ONE LOOP CALCULATION OF THE NN T-MATRIX 


In this section we present the analytical results for the NN T-matrix in one-loop ap¬ 
proximation. We evaluate the contributions at second and third order of the expansion in 
powers of small external momenta (here p, q and m n ). These contributions are divided into 
three classes: vertex and propagator corrections to one-pion exchange, irreducible two- 
pion exchange and iterated one-pion exchange. The latter gives rise to a non-vanishing 
imaginary part in the amplitudes and thus restores unitarity (perturbatively). The possi¬ 
ble local NN-contact terms [||] are not considered here, since to the order we are working 
the corresponding polynomial amplitudes do not contribute to the phase-shifts with L > 2 
and mixing angles with J > 2. 


A. ONE-PION EXCHANGE 

Fig.l shows the one-pion exchange diagram together with (some) one-loop vertex 
and propagator corrections to it. We have evaluated all possible diagrams of this sort 
and found that such one-loop corrections only contribute to mass and coupling constant 
renormalization. No pion-nucleon ’’form factor” is generated by any of these diagrams, 
neither at second nor at third order in the small momentum expansion. Consequently, the 
one-pion exchange diagrams of Fig.l lead just to the following well-known contribution 
to the NN T-matrix, 


W T = 


9^n 


AM 2 (ml + q 2 ) 


( 11 ) 


with g nN , M, m n the physical values of the pion-nucleon coupling constant, nucleon mass 
and pion mass. We also note that in the CM frame there are no 1/M-corrections at any 
order to eq.(ll). It is an exact result of the fully relativistic pseudovector (or pseudoscalar) 
7 rlVlV-vertex evaluated in the CM frame. Finally, we remark that the pionic counter terms 
of 4 4 J in ref. [|T8| and the on-shell (single) nucleon counter terms of C^n i n re f- |[9| do 


not modify the point-like one-pion exchange eq.(ll). 
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Fig.l: One-pion exchange, tree-level diagram and diagrams with one-loop corrections 


B. IRREDUCIBLE TWO-PION EXCHANGE 


The two-pion exchange diagrams are shown in Fig.2. Their contributions to the NN T- 
rnatrix at second order in small momenta come from insertions of vertices and propagators 
of the leading order 7riV-Lagrangian C]/ N alone. All these diagrams can be evaluated in a 
straightforward manner using the heavy baryon formalism of ref. |12[, with the exception 
of the last one, the planar box graph. Here the problem of the so-called pinch singularity 
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arises In the heavy mass limit the two nucleon propagators lead to an expression, 
[(/o + w?)(/o — *h)] _1 with infinitesimal p > 0. The integration contour along the real /o-axis 
is pinched between two nearby poles at ±ir). A loop integral involving this product of 
propagators does not exist since the two poles at infinitesimal distance 2 r) let the result 
diverge as 1/?/. In ref. [|7| is was therefore suggested that one should give up covariant 
perturbation theory and switch over to ’’old fashioned” time-ordered perturbation theory 
in order to avoid nearly vanishing energy denominators. This strategy was then carried 
out in the work of Ordonez, Ray and van Kolck ||. 



Fig. 2: Two-pion exchange 
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diagrams at one-loop order 


The abovementioned problem is of purely kinematical origin. It comes from taking 
the infinite nucleon mass limit M —> oo from the very beginning. In a relativistic cal¬ 
culation with finite nucleon mass M the planar box diagram is well defined, modulo 
ultra-violet divergences which can be handled by dimensional regularization. The planar 
box diagram includes the iterated one-pion exchange which is usually generated by the 
non-relativistic Lippmann-Schwinger equation in the form [|1||]. The non- 

relativistic two-nucleon propagator Q NN is the inverse kinetic energy difference (between 
initial and intermediate state) and thus proportional to the nucleon mass M. It is exactly 
this factor of M in the iterated one-pion exchange which makes the heavy baryon limit 
of the planar box graph divergent and ill-defined. 

The task is to separate off the iterated one-pion exchange proportional to M and 
to find the proper irreducible part which then exists in the infinite nucleon mass limit, 
M —> oo. For that purpose one considers the product of the four relativistic propagators 
entering the planar box diagram and performs the /o-integration via contour methods 
before doing the 1/M-expansion, 


din 


2 vr [(pi - If 


M 2 + irf\[(p 2 + l) 2 - 


M 

[p 2 — (p — l) 2 + irf\ u 2 u >2 


AM 2 i 

M 2 + irj\ [(l + q) 2 — m 2 + irj\ [l 2 
+ + + 0(M-') . 

2(jJ-^CO 2 1 ^2) 


ml + irj\ 


( 12 ) 


Here p± = (E,p), = (-F, —p) and Q 11 = (0, q) are the initial four-momenta of the 

nucleons (see Fig.2) and the four-momentum transfer in the CM frame, uq = \Jl 2 + ml 

and ui 2 = \J\l + q) 2 + ml are the on-shell energies of the two exchanged pions. The factor 
AM 2 was put into the numerator in order to have the correct heavy baryon limit. In the 
hrplane the integrand in eq.(12) has eight simple poles, four below and four above the real 
axis. Two of these poles (one above and one below the real axis) move towards each other 
for increasing nucleon mass M. The integral in eq.(12) can be easily evaluated by closing 
the /o-contour e.g. in the lower half-plane and using residue calculus. The first term in the 
second line of eq.(12), the ’’iterated one-pion exchange”, comes just from that pole which 

















is moving towards another one. The second term, the ’’irreducible two-pion exchange”, 
arises from the two poles of the pion propagators. The fourth pole leads to a strongly 
suppressed contribution of order M~ 5 . The separation into ’’iterated one-pion exchange” 
and ’’irreducible two-pion exchange” as obtained in eq.(12) is not the standard one of 
time-ordered perturbation theory, which gives [2o;f (uq + uo 2)] -1 for the irreducible 
part. Nevertheless, eq.(12) is the correct result and it perfectly agrees with the one of 
time-ordered perturbation theory when performing the systematic 1/M-expansion. Let 
ST = [(p — l) 2 — p 2 )]/2M denote the (single) nucleon kinetic energy difference between 
intermediate and initial (or final) state. From the two irreducible and four reducible time- 
orderings of the planar box diagram one obtains the following expression for the energy 
denominators, 


4u>iU>2 


(uj i + U 2 ) (uJi + ST) (u>2 + ST) 2ST(uj 1 + ST)( u7 2 + ST) 


1 , + ^ 1^2 + UJ 2 - 

2ST uj\uj\ 2ufu)2(uJi + cu 2 ) 


(13) 


including the phase space density factor (4uqc<; 2 ) _1 . 

One can apply the previous method to the crossed box diagram as well and finds as 
a result (in the M —> 00 limit) the same irreducible part as in eq.(12) but with opposite 
sign. In the previous analysis the momentum dependence of the 7rATV-coupling was not 
included. However, in a leading non-relativistic approximation it is of the form a ■ l and 
does therefore not interfere with the Iq- integration. 

As a result one finds that the isoscalar central amplitude Vq adds up to zero for the 
sum of planar and crossed box graphs. For the isovector central amplitude Wc there is a 
further relative minus sign from the isospin factors, which are 3 — 2t\ ■ r 2 for the planar 
box and 3 + 2 ti • r 2 for the crossed box. To the amplitude Wc the two diagrams contribute 
with equal sign. In the spin-spin and tensor channels the situation concerning relative 
signs is just reverse, due to a further relative sign between planar and crossed box graph 
coming from the different ordering of cr-matrices. Altogether, one finds from the graphs 
in Fig.2 the following (renormalized) contributions to the NN T-matrix from irreducible 
two-pion exchange at second order in small momenta, 


W c = 


384t r 2 /. 4 


6ml(15g\ - 6 g\ - 1) + g 2 (23^ - 10^ - 1) 


m n 

ln - 


+4m?(4 g\ + + 1) + -jr(5g\ + 2 (ig 2 A + 5) 

o 


+ 


4m 2 (5 <7 — 4 g\ — 1) + q~(23g\ — 10 g\ — 1) + (L I L(q) \ , 


4m^ + q 2 \ 


V T = = 


34 

647 f 2 f‘ 




(14) 

(15) 


with the logarithmic loop function 


L(q) 


w w + q 

— In-, 

q 2 m n 



(16) 


The NN-amplitudes not specihed in eqs. (14,15) are zero. The first two lines in eq.(14) for 
Wc and the constant terms in eq.(15) for Vt do not contribute to phase shifts with L >2 
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and mixing angles with J > 2. Such polynomials in q 2 drop out of the projection formulas 
eqs.(6-9) by orthogonality of the Pj(z). Therefore the value of the (a priori arbitrary) 
renormalization scale A is irrelevant in this work. 

Next we consider the two-pion exchange contributions at third order in the small 
momentum expansion. These arise from the one-loop diagrams in Fig.2 with exactly 
one insertion from the next-to-leading order 7TTV-Lagrangian The contributions are 

either proportional to the low-energy constants c, or they carry a suppression factor 1/M. 
Again, all but the planar box graph can be evaluated in straightforward manner using the 
heavy baryon formalism 0 - The planar box graph requires a special treatment in order 
to obtain correctly all 1/M-corrections. First, the expansion of the propagator integral, 
eq.(12), in powers of 1/M has to be performed one order further. This gives rise to new 
irreducible pieces (proportional to 1/M) and a relativistic correction factor, l — p 2 /2M 2 ~ 
M/FT, to the iterated one-pion exchange. Secondly, one has to expand in powers of 1/M 
the product of pseudovector 7rTVTV-vertices and (J> + M)/2M-factors from the nucleon 
propagators sandwiched between in- and out-going Dirac spinors. At order 1/M, only 
terms proportional to Z 0 and Zq are generated this way. The integral analogous to eq.(12) 
with an additional factor Zo or Zq in the numerator gives zero, 0 + 0(M -2 ). At order 1/M 2 
terms independent of Z 0 but proportional to Z • (2 p — l ) appear. Together with eq.(12) 
these give rise to a further 1/M-correction to the irreducible two-pion-exchange. With 
these techniques and the rules of dimensional regularization to evaluate the remaining d 3 l- 
integral one obtains all 1/M-corrections coming from the planar box graph. As a check 
we applied the same technique to the crossed box and triangle graphs. In all cases we 
found perfect agreement with the results obtained independently using the heavy baryon 
formalism. This shows that our method of doing first the Z 0 -integral and then expanding 
in 1/M, which is necessary for a systematic evaluation of the planar box diagram, is 
indeed correct. Putting all pieces together, we finally arrive after tedious calculation at 
the following irreducible two-pion exchange contributions to the NN T-matrix at third 
order in small momenta, 




3 g\ 


{ 4 (ci - cz)ml - + 3 q 2 ) 
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(17) 


(18) 
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( 21 ) 


Vso 

kE so 
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+ 4 ^) (4m ' + ~ + 3q ^ ’ 

^ j I {m n + (2ml + q 2 )A(q )} , 

{ra*- + (4 m 2 + g 2 )A(g)} . 


QAttM f. 

A( i - <fl) 


( 22 ) 

(23) 


Note that the Ci-part of W in eq.(17) is proportional to the one-loop contribution to 
the nucleon scalar form factor, namely — 8 ci /“ 2 cqv(—g 2 )/ oop j^ 0 | and the low-energy con¬ 
stant C 3 is related to the so-called nucleon axial polarizability |[T 2 | . |TT[| , a A = - 2 c 3 /fl - 
g\rn n (4:8g\ + 77)/(3847 rff) = 6.44 fm 3 which is governed by virtual A-isobar excita¬ 
tion. Similarly, the C 4 -pa.rt of Wt in eq.(21) is proportional to the one-loop contribu¬ 
tion to the nucleon isovector magnetic form factor (normalized to 1 + k v — K n ), namely 
— {ci/AM fl) Gm(— q 2 )i oop - The low-energy constant C 2 is absent in the above expressions. 
In the actual calculation one sees that this comes from the fact that the energy transfer in 
the CM frame is zero, g 0 — 0. We note that all contributions at third order in small mo¬ 
menta are finite (free of divergences). This is consistent with the fact that no local counter 
terms (analytic in the quark mass) exist at this order of the small momentum expansion 
to cancel divergences. The spin-orbit amplitudes Vso , Wso generated by chiral two-pion 
exchange are proportional to 1/M. This is different from the OBE-model, where scalar 
and vector meson exchange lead to spin-orbit terms proportional to 1/M 2 . Furthermore, 
to the order we are working, there is no quadratic spin-orbit contribution from irreducible 
graphs. Such terms are suppressed by a factor 1/M 2 . Already a short glance at the ana¬ 
lytical formulas, eqs. (14-23) for the two-pion exchange contributions shows that these are 
all increasing (in magnitude) with increasing momentum transfer q. This is in contrast to 
the OBE-model, where by construction all NN-potentials decrease as a function of q. The 
growth in q of the chiral two-pion exchange potentials is a consequence of the Goldstone 
boson nature of the pion. A Goldstone boson has (and induces) weak interaction at low 
energies and stronger interaction at higher energies (see also the explicit expressions for 
the one-loop 7 rlV-scattering amplitude in ref. [jTIJ and the one-loop 7 r/r-scattering ampli¬ 
tude in ref. [|ll|, which share the same property). However, in many model calculations, 
this feature is altered by the introduction of phenomenological cut-offs, form factors and 
the like. 


C. ITERATED ONE-PION EXCHANGE 

The planar box graph includes the iterated one-pion exchange, which is enhanced 
by a factor M. Even though we counted the diagram naively as second order in small 
momenta, the actual contribution to the NN T-matrix is of first order in small momenta. 
The occurrence of such a large scale enhancement factor M in iterated diagrams changes 
the correspondence between the loop expansion and the small momentum expansion for 
processes involving two (or more) nucleons. Ladder diagrams with n loops contribute 
already at n-th order in the small momentum expansion, whereas the naive counting 
rules would classify them as order 2 n. The iterated one-pion exchange including the 
relativistic correction factor M /E has the following integral representation, 
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The isospin factor 81 — 9 suppresses the contribution to (/ = Instates by a factor 1/9 
compared to (/ = 0)-states. The resulting three-dimensional integrals can be evaluated 
in closed form, and the contributions to the NN T-matrix are expressed in terms of the 
following set of complex-valued functions, 
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(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


with w and A(q) given in eqs.(16,18). Note that the functions Gg 2 , 3 ( 1 °, q) are not singular 
at z = — 1 or q 2 = 4p 2 . In terms of these complex functions the contributions to the NN 
T-matrix read, 
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= 3 ^^ 4 { g, o(p, g) + 4Gi(p, g) + 4G 3 (p, g)} . (34) 


3 ^ 64vr Efl 

4 /\/r2 


flere, for the first time a non-vanishing quadratic spin-orbit contribution is found. Using 
dimensional regularization, we find that the iterated one-pion exchange is finite and has 
no divergences. Note that the contributions of iterated one-pion exchange, eqs. (31-34), 
depend on both kinematical variables p and g, whereas irreducible two-pion exchange, 
eqs. (14-23), depends only on the momentum transfer q. To the order we are working 
the iterated one-pion exchange is the only contribution to the NN T-matrix with a non¬ 
vanishing imaginary part and thus it restores unitarity to this order. In partial waves this 
imaginary part is given by the square of the phase-shift calculated from one-pion exchange, 
eq.(ll); in the case of the mixing triplet states with L — J ± 1 it is actually a quadratic 
expression in these phase shifts and the mixing angle. Therefore the imaginary part of 
iterated one-pion exchange allows for a numerical check of unitarity, which actually serves 
as a check of both, the analytical expressions for iterated one-pion exchange, eqs. (31-34) 
and the projection formulas, eqs. (6-9). We have found that in all cases this numerical 
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check of unitarity works at an accuracy 10 -4 or better. Furthermore, we have verified 
analytically that the optical theorem is satisfied. It equates the imaginary part of the NN 
T-matrix in forward direction (only Vq and Wq in eq.(31) contribute) to pE/M 2 times 
the (unpolarized) total NN-cross section calculated in (non-relativistic) one-pion exchange 
approximation without crossed nucleon lines. 

Since each iteration of one-pion exchange brings a further enhancement factor M, up 
to three iterations will contribute to the NN T-matrix calculated completely to third order 
in small momenta. This corresponds to two- and three-loop diagrams which are extremely 
difficult to evaluate. Numerically, we find that in most cases iterated one-pion exchange 
is small compared to irreducible two-pion exchange. We can therefore assume that the 
higher order iterations can be neglected. 

At zeroth order in the small momentum expansion, not only one-pion exchange eq.(ll) 
contributes to the NN T-matrix, but also contact terms of the (Fierz-antisymmetric) form, 

3 B + B' Ti ■ r 2 — (2 B + B') d\ • cr 2 — £> fj • r 2 d\ ■ cf 2 . (35) 

In tree approximation eq.(35) gives a contribution to the two S-wave amplitudes only. At 
one-loop order the iteration of this contact interaction and diagrams with an additional 
pion exchange occur. We have checked that the one-loop diagrams involving the lowest 
order contact interaction eq.(35) do not contribute to the phase-shifts with L > 2 and 
the mixing angles with J > 2. Therefore, we do not need to specify the constants B, B' 
in eq.(35). 

This completes the discussion of the NN T-matrix in one-loop approximation. We 
fold that if we restrict ourselves to the phase-shifts with L > 2 and mixing angles with 
J > 2 (peripheral partial waves), not even a single adjustable parameter shows up. To the 
order we are working this peripheral part of the NN T-matrix involves only well-known 
physical parameters such as g-^N, m n) M and three low-energy constants 01 , 3,4 which 
are well determined from 7 rWscattering data. Such features allow for a clean test of chiral 
symmetry in NN-interaction. 


V. RESULTS AND DISCUSSION 


I 11 this section we present and discuss our results for the phase-shifts with L > 2 
and mixing angles with J > 2 up to nucleon laboratory kinetic energies of T; a b = 280 
MeV, the AW 7 r°-threshold. At such an energy, momentum transfers up to q = 5.25 m n = 
725 MeV are involved, which are quite large for the application of chiral perturbation 
theory. Despite this fact, we present here results up to the AW 7 r°-threshold, in order 
to demonstrate where our model independent predictions are in agreement with existing 
data and where deviations show up. For the sake of completeness, we summarize again 
the values of the physical parameters used in the calculation: = 92.4 MeV, g n N = 13.4 

as determined from 7 rV-dispersion analysis |L3|, m n = 138 MeV, M = 939 MeV for 


the average pion and nucleon mass and the central values of the low-energy constants, 
Ci = —0.9 GeV -1 , c 3 = —5.3 GeV -1 and C 4 = 3.6 GeV ^ 1 


m 
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A. D-WAVES 


The D-wave phase shifts and mixing angle e 2 are shown in Fig.3. The dashed curve 
corresponds to the one-pion exchange approximation and the full curve includes in addi¬ 
tion two-pion exchanges (irreducible two-pion exchange and iterated one-pion exchange). 
The dotted and dashed-dotted curves represent the empirical NN phase shift analyses of 
ref. and ref. [21] (as far as available), respectively. In all cases the model independent 


two-pion exchange corrections go into the proper direction, but they tend to be too large, 
already at small Ti ab . The one-loop predictions for the D-wave phase shifts deviate appre¬ 
ciably from the data. An exception is the mixing angle 62 which is well reproduced up to 
Ti a b = 120 MeV. The 27r-exchange tensor force in this channel has the correct sign. How¬ 
ever, it grows too strongly with energy and even overcompensates the l7r-exchange tensor 
force at the AUV 7 r 0 -threshold. The 3 D\ phase shift is also in fair agreement with the data 
up to Tiab = 200 MeV. This results from an almost complete cancelation of irreducible 
two-pion exchange and iterated one-pion exchange contributions. Above Ti ab = 200 MeV 
this cancelation mechanism does not work anymore. The good agreement between the 
data and the one-pion exchange approximation in the 3 D 1 partial wave is contrary to 
expectations and appears accidental, since in G-waves (see later) there are still sizeable 
differences between data and one-pion exchange approximation. From all this we conclude 
that one-pion and two-pion exchange alone are not sufficient to describe the dynamics in 
the NN D-waves [|5j (see also the discussion of the coordinate potentials in the next sec¬ 
tion). Obviously, strong cancelations between 27r-exchange and 37 r-exchange together with 
shorter range effects are taking place here. The latter ingredients of the NN T-matrix go 
beyond the one-loop approximation. 


B. F-WAVES 


The F-wave phase shifts and the mixing angle 63 are shown in Fig.4. Again, the 
27r-exchange corrections to l 7 r-exchange go into the right direction. The predictions for 
the phase-shifts in the 1 F 3 and 3 F 3 partial waves are in good agreement with existing 
data up to Ti ab = 180 MeV. At larger energies the 27r-exchange effects become again 
too strong. In the 3 F 2 partial wave we find a correction to l 7 r-exchange opposite to the 
trend of the data of ref. |2l|. The latest fit of the 3 F 2 phase shift from the VPI group 


Q (not shown) is very similar to the one of ref. The mixing angle e 3 , however, is in 


perfect agreement with the data of ref. [|4|,^1[] for all energies up to the VAUffi-threshold 
at Ti a b = 280 MeV. To this observable irreducible 27r-exchange and iterated l 7 r-exchange 
contribute with almost equal strength. Note that the empirical F-wave phase-shifts are 
quite small, namely less than 4.5°. It would be very interesting to see how the present 
discrepancies between the chiral NN F-wave phase shifts and the existing analyses reflect 
themselves in a direct comparison with scattering data, i.e. differential cross sections and 
polarization observables. Unfortunately, we can not present such a direct comparison 
here. 
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C. G-WAVES 


The G-wave phase shifts and the mixing angle e 4 are shown in Fig.5. The chiral 
predictions are in good agreement with data for all four partial wave phase shifts and 
the mixing angle up to Ti ab = 220 MeV. Note that the differences between l 7 r-exchange 
and data are still sizeable in the 1 Gb and 3 Gs partial waves. The model independent 
chiral 27r-exchange corrections have indeed the correct sign and magnitude to close this 
gap between data and l 7 r-exchange approximation. I 11 the 1 G± partial wave the correction 
comes mainly from irreducible 27r-exchange, whereas in the 3 G§ partial wave (with 1 = 0) 
both irreducible 27r-exchange and iterated l 7 r-exchange contribute with roughly equal 
strength. For the 3 G 3 and 3 G 4 phase shifts and 64 the 27r-exchange corrections are quite 
small. Nevertheless, these small effects improve the agreement between data and chiral 
predictions. 


D. H-WAVES 

The H-wave phase shifts and the mixing angle 65 are shown in Fig. 6 . The corrections 
due to 27r-exchange are already quite small and they follow the trend of the data. An 
exception is the 3 H 4 phase shift where we find a very small positive contribution. However, 
the empirical 3 7 / 4 phase shifts are very small, less than 0.5°, and the difference between 
the existing phase shift analyses and the chiral prediction may be insignificant if one 
compares directly the NN-scattering data. 


E. I-WAVES 


The I-wave phase shifts and the mixing angle 6$ are shown in Fig. 7. Again, the 
corrections coming from 27r-exchange are quite small and they follow the trend of the 
data. Altogether, we find a faster convergence to the l 7 r-exchange approximation in the 
high angular momentum partial waves as obtained in the existing NN phase-shift analyses. 
The calculation in chiral perturbation theory presented here is most reliable in the high 
angular momentum partial waves. These are the ones which probe the long and medium 
range parts of the NN-interaction which are built in here in a model independent way via 
17T— and chiral 27T-exchange. In the existing phase shift analysis the (small) effects in high 
angular momentum partial waves may result mainly from effective potentials [j|,^l|] fitted 
to lower partial waves with only limited direct constraints from the NN-scattering data. 
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FIGURES 




Fig. 3: D-wave NN phase shifts and mixing angle 62 versus the nucleon laboratory ki¬ 
netic energy Ti a b. The dashed curve corresponds to the lTT-exchange approximation and full 
curve includes chiral 2n-exchange as well. The dotted and dashed-dotted curves represent 
the empirical NN phase shift analyses of ref. a and ref. respectively. 
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Fig.4: F-wave NN phase shifts and mixing angle 63 versus the nucleon laboratory 
kinetic energy Ti ab . For notation see Fig. 3. 
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Fig. 5: G-wave NN phase shifts and mixing angle 64 versus the nucleon laboratory 
kinetic energy T tab . For notation see Fig. 3. 
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Fig. 1: I-wave NN phase shifts and mixing angle e§ versus the nucleon laboratory kinetic 
energy Ti ab . For notation see Fig. 3. 
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VI. COORDINATE SPACE REPRESENTATIONS 


In order to interpret the effects of irreducible chiral two-pion exchange and to compare 
with phenomenological approaches it is desirable to have a coordinate space representation 
of the corresponding potentials. However an ordinary inverse Fourier-transform of the 
momentum space functions given in section IIB is not possible, because of their growth 
with q. In this case the coordinate space representation of the finite range parts can be 
obtained in the form of a continuous superposition of Yukawa functions [jl|. The mass 
spectrum entering in this representation is given by the imaginary part of the momentum 
space functions analytically continued to time-like momentum transfer, q = 0 + — i p. For 
the isoscalar central potential in r-space one has then the following spectral representation, 

Vc{r) = -7TT [ dppe~^ r lmV c (-ip) . (36) 

zn-r J2m.Tr 


Using Im A(—ip) = (n/4:p)9(p — 2m n ) for the relevant loop function (note that due to the 
coalescence of normal threshold /i 0 = 2?7and anomalous threshold p c = J 4 m 2 


mf/M 2 

= 1.995771^ in the heavy nucleon limit, this imaginary part starts with a non-zero value 
at threshold) one finds the following closed expression for the isoscalar central potential, 


V c {r) = 


3 g\ 


t-2x 


327t 2 /4 r 6 


(2ci + 


M 

16 M 


x 2 (l + x ) 2 + 


9a % 5 

32 M 


+ (c 3 + y^) (6 + 12x + 10x 2 + 4x 3 + X 4 ) | 


(37) 


with the abbreviation x = m n r. Since ci^ + 3g\/16M < 0, this potential is attractive 
for all distances r as expected. Its asymptotic behaviour for large r is a modified Yukawa 
function with twice the pion mass, e~ 2m ' nr / r 2 , neglecting the small 1/M-terms. At a 
distance r = 1 /m n = 1.43 fm one finds an attraction of Vc(1.43fm) = —35.9 MeV. This 
number is twice as large as the attraction produced by the fictitious cr(550)-boson of the 
Bonn potential [|I].[5| . 




Fig. 8: The isoscalar central potential Vc(r ) generated by irreducible chiral two-pion ex¬ 
change in r-space (full line). The dashed lines show the phenomenological a + oj and 
a-contributions. 


In Fig.8 we show the isoscalar central potential Vc(r ) for 1 fm < r < 2 fm and 2 fm 
< r < 4 fm. Our results are qualitatively in agreement with the recent work of ref. |23| 
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but their isoscalar central NN-potential is even somewhat more attractive. The terms 
proportional to in eq.(37) agree with the earlier work of [^4|] up to the isospin factor 
3 which was inadvertently omitted in |24[]. Furthermore it is interesting to observe that 
in the chiral limit m n = 0 the isoscalar central two-pion exchange potential shows a 1 /r 6 
fall-off which is typical for a non-relativistic Van-der-Waals potential. 

The spectral representation of the isovector tensor potential (accompanying the stan¬ 
dard tensor operator 3<fi • f dy ■ r — a\ ■ oq) takes a similar form, 


W r (r) = 


d/ilie Mr (3 + 3 fir + p 2 r 2 ) ImkF r (—i/i) . 


67 r 2 r 3 J 2 m, 

The two-pion exchange isovector tensor potential reads explicitly, 


W^\r) = 


9a e ~ 2x f / 1 


48^2/4 r « {-F + m )(l + x)(3 + 3x + A) 

(36 + 72x + 52x 2 + 17x 3 + 2x 4 )| . 


9a 

32 M' 


(38) 


(39) 


Since C 4 > 0 it is attractive and therefore cuts down the repulsive isovector tensor potential 
due to one-pion exchange 


W^\r) 


9nN e x 
48itM 2 r 3 


(3 + 3a; + x 2 ) . 


(40) 


At a distance r = 1 /m n = 1.43 fin the two-pion exchange isovector tensor potential is 
—22.5% of the one-pion exchange isovector tensor potential. Such a reduction is indeed 
required by phenomenology |I|. In Fig.9 we show the reduced isovector tensor potentials 
(m 7 r r) 3 PFr(r) for distances 1 fin < r < 4 fin. Note that already at r = 2.5 fm the two-pion 
exchange effects start to become visible. The 1/M-terms in eq.(39) are very small. 



Fig.9: The reduced isovector tensor potentials (m n r) 3 Wt{t) in r-space. 

Furthermore, we give the two-pion exchange isoscalar tensor and isovector central 
potentials in coordinate space. Using Im L(—ifj) = —(ti/2/j,)^/j, 2 — 4m 2 for the relevant 
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loop function, one finds after some calculation that these can be expressed in terms of 
two modified Bessel-functions, 


V T {r) = 


g\ m n 


1287 r 3 /^r 4 
3Tre~ 2x 


— 12xK 0 (2x) — (15 + 4a; 2 ) R\(2x) 
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W c {r) = 


8 Mr 
rrw 


1287 r 3 /^r 4 


(12a; 1 + 24 + 20a; + 9a; 2 + 2a; 3 ) | , 

[1 + 24(5 + 2x 2 ) - 4(23 + 12a; 2 )] K x (2x) 


(41) 


+a;[l + 104 — 4(23 + 4a; 2 )] Kq{2x) + 


2 —2r 

9A* e 
4: Mr 


2(3 g\ - 2) 


•(6a; + 12 + 10a; + 4a; + x 6 ) + g A x{2 + 4a; + 2x + 3a; 3 ) 


(42) 


The asymptotic behaviour for large distances r is e~ 2mirr /r 5 / 2 for Vr(' r ) and e _2m7rr /r 3 / 2 
for Wc(r), neglecting the 1/M-terms. This long range behaviour is determined by the 
box graphs (proportional to g\). In the chiral limit m n = 0 both display a 1/r 5 powerlike 
fall-off. Since the modified Bessel-functions K 0t i(2m n r) are positive and g A > 1 one sees 
that Vt( r) and Wc(r) are both negative throughout. The finite range part of the two- 
pion exchange isovector central potential is attractive and not repulsive as one would 
expect from analogies with the p-meson. It is very interesting to see how this unexpected 
attraction comes about. In the bubble and triangle diagrams (proportional to g A and g A ) 
of Fig.2 a two-pion state with the quantum numbers of the p-meson is emitted from one 
nucleon. As can be seen from the explicit expression in eq.(42) these processes indeed 
lead to an isovector central repulsion. However, the box diagrams (proportional to g A ) 
are much more attractive and overcompensate the expected isovector central repulsion. 
The isovector central and isoscalar tensor potentials are shown in Fig.10. We remark that 
for these the 1/M-terms are not negligible. 



Fig. 10: The isovector central potential Wc(r) and the isoscalar tensor potential Vr(r) 
generated by irreducible chiral two-pion exchange in r-space. 
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At r = 1 fm they lead to an almost 50% reduction of the leading term, while at larger 
distances this reduction gets smaller. 

Two-pion exchange spin-spin potentials in coordinate space (accompanying the oper¬ 
ator ai ■ a 2 ) have an analogous spectral representation 

~ 1 r°° 

v s{r) = -— r dfj, n 3 e~^ r lmVr(-ifi) . (43) 

37r z r J 2 m- n 


We find from the chiral two-pion exchange the following isoscalar and isovector spin-spin 
potentials, 


Y (y) — 9A m Tr J 1 /O I Q.v.S'i 
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Both are positive and thus add to the one-pion exchange spin-spin potential 


(44) 


(45) 


W^(r) - g " Nm * 6 X 
Ws {r) ~ 487rM 2 r ' 


(46) 


In Fig. 11 the spin-spin potentials due to one-pion-exchange (isovector) and two-pion ex¬ 
change (isoscalar and isovector) are shown for 1 fin < r < 2 fin. The isoscalar spin-spin 
potential is again reduced appreciably by the 1/M-terms in eq.(44). We also note that the 
terms proportional to g\ involving the modified Bessel-functions Kq^{2x) in eqs. (41,42,44) 
agree with the earlier calculation of ref. 
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r [fm] 


Fig. 11: Spin-spin NN-potentials due to one-pion exchange (dashed line) and two-pion 
exchange (full lines) in r-space. 
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Spin-orbit potentials in coordinate space (accompanying the standard spin-orbit op¬ 
erator — |(<7i + cr 2 ) ■ (r x V r )) have a spectral representation, 


Vso(r) = -^-g / d/i /i (1 + /ir) e Im V so (-in ) . (47) 

7T r j ‘Im-r: 

We hnd from chiral two-pion exchange the following isoscalar and isovector spin-orbit 
potentials, 


Kso(r) 

^so(r) 


34 e- 2a: 

647t 2 Mr 6 

4 M ~!) e ~ 2 * 

327 r 2 Mff r 6 


■(1 + x) (2 + 2x + a; 

(1 + x ) 2 . 


(48) 

(49) 


At r = 1/mjr = 1.43 fin these finite range two-pion exchange spin-orbit potentials are 
Vso(l/m n ) = —1.97 MeV and Wsoi)-h^n) = 0.22 MeV. The isovector spin-orbit potential 
hhsoM is very small and repulsive. However, the isoscalar spin-orbit potential Vso{ r ) 
from irreducible two-pion exchange is quite large and attractive and it turns out to be in 
very good agreement with the sum of the phenomenological c(550)- and o;(782)-exchange 
contributions. In Fig. 12, these spin-orbit potentials are shown for distances 1 fin < r < 2 
fin. Finally, we note that the expressions for the two-pion exchange spin-orbit potentials 
eq. (48,49) agree with the earlier calculation of ref. |2^|. Actually, we agree on all terms 
with earlier calculations, except those proportional to g\/M coming from the planar 
box graph and contributing to the central, spin-spin and tensor potentials. This type 
of relativistic correction is only well defined if the relativistic correction to the iterated 
one-pion exchange is specified (as done here in eq.(24)). 



Fig. 12: Isoscalar and isovector spin-orbit NN-potentials of two-pion exchange in r-space 
(full curves). The dashed line shows the isoscalar spin-orbit potential from a- and u- 
exchange Ml- 

As a summary we hnd, that all irreducible two-pion exchange NN-potentials are of 
van-der-Waals type. They show an asymptotic exponential behaviour, e~ 2rrinr /r n , n = 
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|,2,|,3, 4, with a decay length (2m 7r )~ 1 ~ 0.7 fm. Near the origin r = 0 they become 
singular as r ~ 5 or r -6 . The strongly attractive parts of the potentials Vo, Id's and Wt 
below r = 1.5 fm are responsible for the overestimate of the D-wave phase shifts and 62 in 
Fig.3. Wave-functions of D-states still probe to an appreciable amount the NN-potentials 
below r = 1.5 fm, which in our calculation come out too strong in most cases. In order to 
improve the D-wave observables one has to add either repulsive zero-range counter terms 
or an explicit (782)-meson. 

Finally, we have to spell out a warning in order to avoid confusion. The coordinate 
space potentials calculated here and the momentum space functions of section IIB are 
not related through each other by ordinary Fourier- and inverse Fourier-transformation. 
Both of them are too singular at r = 0 or q —> 00 that these transforms would exist. 
The polynomial pieces (a + (3q 2 ) of the momentum space functions generate potentials 
proportional to zero-range delta-functions and Laplacians of delta-functions with diver¬ 
gent constants. On the other hand the Fourier-transform of a potential as singular as 
r -5 Qr r -6 near the origin has a divergent polynomial piece in it. Only in a loose sense, 
i.e. modulo zero-range <5 3 (r )-terms and linear polynomials in q 2 , are the coordinate space 
potentials (presented in this section) and the momentum space functions of irreducible 
two-pion exchange (section IIB) related by Fourier-transformation. For the observables 
we have considered in this work, the phase shifts with L > 2 and mixing angles with 
J > 2, these polynomial pieces and delta-function terms do not play a role at all, as we 
have mentioned repeatedly. 


VII. COMPARISON WITH ONE-BOSON-EXCHANGE 


A detailed comparison between one-boson exchange and (irreducible) chiral two-pion 
exchange in momentum space is not really meaningful because of the different dependences 
on the momentum transfer q. Nevertheless, we find it instructive to compare at least the 
overall strengths of both, i.e. the values at q — 0, denoted by an overbar. These values at 
q = 0 include also zero-range effects generated by irreducible chiral two-pion exchange. 

For the isoscalar central part one finds from irreducible chiral two-pion exchange and 
<7 + tu-exchange at q — 0, 


T7 (27r) 
v c 


tjD) , tjG) 


V + V 


c 


f 6 9a \ _ 

16tt/ 4 V 5 Cl C3 64M/ 

A _ A = (295 - 217) GeV 


= 79 GeV“ 2 , 

" 2 = 78 GeV" 2 , 


(50) 


using g 2 /47r = 7.1 and g 2 /47 t = 10.6 [1],[5|]. The isoscalar central attraction produced 
by chiral two-pion exchange is stunningly close to the sum of the individually large phe¬ 
nomenological cr-attraction and cu-repulsion. Note that the isoscalar central attraction due 
two-pion exchange is a chiral symmetry breaking effect ~ m 3 . The following interpreta¬ 
tion of the numerical result, eq.(50), is tempting. In the chiral limit the phenomenological 
a-attraction and cu-repulsion will cancel each other entirely. On the other hand, it has 
been shown recently |22J that in this case, Vo\r) + Vo"\r) ~ 0, the so-called pseudospin 
symmetry observed in nuclear shell model spectra, can be explained very naturally. This 
suggests a strong connection between chiral symmetry of QCD and pseudo-spin symme¬ 
try in nuclei. Of course, more detailed and quantitative work is necessary to put this 
conjecture on a firm basis. 


26 





We have also compared irreducible chiral two-pion exchange with OBE in other chan¬ 
nels. No clear relation between OBE (a, p, uj) and chiral two-pion exchange can be deduced 
from this comparison. The physics behind both is also rather different. Whereas the OBE 
parametrizes effects of (strongly) correlated and resonant multi-pion exchange, the chiral 
two-pion exchange (to one-loop) involves only soft pions without self-interaction. 


VIII. SUMMARY 

In this work we have gone one step further in a model independent description of 
the low-energy NN-interaction. In addition to the one-pion exchange, which is known to 
dominate the high angular momentum partial waves, we have included two-pion exchange 
processes based on the most general chiral effective Lagrangian. We have performed a one- 
loop calculation of the NN T-matrix using covariant perturbation theory and dimensional 
regularization throughout. We solved the so-called pinch singularity problem occuring in 
the planar box diagram by subtracting the iterated one-pion exchange proportional the 
nucleon mass M. In fact we calculated here all one-loop contributions to the NN T-matrix 
at second and third order in the small momentum expansion. 

No pion-nucleon form factor is generated to modify the point-like one-pion exchange. 
Restricting ourselves to peripheral NN-phases (phase shifts with orbital angular momen¬ 
tum L > 2 and mixing angles with J > 2) the calculation does not involve a single 
adjustable parameter. This feature allows for a clean test of chiral symmetry in the two- 
nucleon system. Comparing to existing NN phase shift analyses we find partial agreement 
in the D-waves up to T) afc = 150 MeV but also appreciable deviations. In particular, we 
End that the isoscalar central attraction generated by chiral two-pion exchange is far too 
strong for distances 1 fm < r < 2 fm. This points towards the importance of shorter 
range effects in the D-waves. 

Increasing the orbital angular momentum L the agreement between the chiral predic¬ 
tion and the data improves gradually also for larger energies up to the Ah/V 7 r°-threshold, 
Ti a h = 280 MeV. There is indeed a kinematical window where chiral symmetry governs 
the nucleon-nucleon interaction. Examples for this are the mixing angle 63 and the 1 G 4 
and 3 G ,5 phase shifts, where the chiral corrections close the sizeable gaps between the 
one-pion exchange approximation and the data. In the high angular momentum partial 
waves (L > 5) we find a faster convergence to the one-pion exchange as obtained in the 
empirical NN phase shift analyses. 

The calculation presented here incorporates via one- and two-pion exchange the long 
and medium long range components of the NN-interaction in a model independent way. 
Therefore we propose the chiral NN phase shifts with L > 3 and mixing angles with 
J > 3 to be used as input in a future NN phase shift analysis. This way one could 
find out whether the small differences to existing phase shift analyses are significant in 
a direct comparison with NN-scattering data, namely angular distributions of differential 
cross sections and polarization observables. 

We have given explicit expressions for the isoscalar/isovector central, spin-spin, tensor 
and spin-orbit potentials in coordinate space generated by irreducible chiral two-pion 
exchange. They are of van-der-Waals type with an exponential asymptotic behaviour, 
e — 2 m w r j r n.' ^he isoscalar central potential is about twice as attractive as the one of the 
(T(550)-boson in the Bonn potential. Contrary to expectations the two-pion isovector 
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central potential is attractive, because the contribution of the box graphs is numerically 
dominant. 

The one-pion isovector tensor potential is reduced by the two-pion exchange contribu¬ 
tion. The sum of both is in agreement with the phenomenological Paris potential above 
r = 1.2 frn. The two-pion exchange isoscalar spin-orbit potential is surprisingly close to 
the sum of the phenomenological (x(550)- and cu(782)-exchange contributions. 

It is needless to say, that the chiral approach to NN-interaction presented here can 
(at the moment) not substitute for the successful phenomenological models as a whole. 
In nuclear physics applications the low partial waves (S, P, D) are the more relevant ones 
and chiral symmetry may be of minor importance for these. A combination of chiral 
dynamics together with short range contributions e.g. from ^-exchange might help and 
will be explored further. 
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APPENDIX: ANTI-SYMMETRIZED NN T-MATRIX 

Even though the (direct) NN T-matrix of the form eq.(4) is sufficient to calculate 
all observables it has one deficit. It does not incorporate the Pauli exclusion principle 
and therefore it has non-vanishing matrix elements in the Pauli-forbidden NN-states with 
even I + S + L. This deficit is cured by subtracting the T-matrix which results from the 
exchange of the two out-going nucleon lines, 


T N n — A[T? 


NN\ 


(51) 


The exchange operation A involves a left-multiplication with the isospin exchange oper¬ 
ator (1 + ri • t 2 )/2 , a left-multiplication with the spin exchange operator (1 + a\ ■ <? 2 )/2 
and a substitution p' —> —p' and z —> — z in the scalar functions Vc(p, z ),..., Wq( p, z ). 
Evidently, twofold application of the exchange operation A is the identity operation, 
A o A = 1, and thus T^n — A\Tnn] is indeed anti-symmetric. The exchange T-matrix 
A[T N n] has the same decomposition with respect to the ten isospin-spin operators as T NN 
in eq.(4), only the scalar functions Vc(p, z ),..., Wq( p, z ) are replaced by different ones, 
called A[Vc(p, z)],.,..., A[Wq(p, z)\. The latter have the following explicit form (suppress¬ 
ing the arguments p and z), 


A[V C ] 

A\W C ] 

A[V S ] 

A\W S ] 

A[V T ] 

A[W t \ 


1 

4 

1 

4 

1 

4 

1 

4 

1 

2 

1 

2 


Pc + SW C + 3P S + 9W S + q\V T + 3 W T ) +p\l- z 2 ) (V Q + 3 W Q ) 
V c - W c + 3W - 3W S + q 2 (V T - W T ) + p\ 1 - ^ 2 )(Pq - W Q j ^ 
V c + ZWc -V s - 3W S + q 2 (V T + 3 W T ) + p\z 2 - 1)(P Q + 3 W Q ) 
Vc -Wc-Vs + W s + q 2 {V T - W T ) + p\z 2 - 1){V Q - W Q ) 1 ^ 


rZ — 1 

■ 2+1 
rZ — 1 
-Z+l 


(V T + 3 W T ) 


— Z 


(' V T - W T ) 


) 


5 



A[V S o] — x ~ Vso — 3 Wso , 

Z J £ 

A[Wso] = ^\Wso-Vso\ 

Z 1 J Z—>~Z 

AVq] = ^[vq + 3W r Q - pi(rTT)^ + . 

■W = \[Vq ~ W Q - -^TJ)(Vt - Wt)]^ , (52) 

with the notation [/(p, z)\ z ->- z = f(p , —z). The matrix elements of A[T NN ] in the LSJ- 
basis are readily evaluated replacing Vc, • • •, Wq by M[V0], • • •, ,4[Wq] in eqs.(5-9). Mak¬ 
ing a further substitution z —> — z and using the property Pj(—z ) = (—1 ) J Pj(z) of the 
Legendre polynomials, one finds that the original matrix elements of Tnn are reproduced 
up to an important sign, 

{L'SJ\A[T nn ]\LSJ) = (-l) /+s+i (L'SJ\T nn \LSJ) . (53) 

The anti-symmetrized T-matrix T NN — A[T NN ] has indeed the desired property, that all 
its matrix elements in the Pauli-forbidden NN-states (with even / + S + L) vanish identi¬ 
cally. Imposing the Pauli principle from the very beginning is important if one considers 
polynomial counter terms in order to avoid redundance in the low-energy constants (see 
e.g. eq.(35) for an S-wave counter term with only two independent low-energy constants 
B, BO- 
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